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Theoretical Design and FPGA-Based 
Implementation of Higher-Dimensional Digital 

Chaotic Systems 

Qianxue Wang, Simin Yu, Chengqing Li, Jinhu Lii, Xiaole Fang, Christophe Guyeux, and Jacques M. Bahi 


Abstract —Traditionally, chaotic systems are built on the do¬ 
main of infinite precision in mathematics. However, the quan¬ 
tization is inevitable for any digital devices, which causes dy¬ 
namical degradation. To cope with this problem, many methods 
were proposed, such as perturbing chaotic states and cascading 
multiple chaotic systems. This paper aims at developing a novel 
methodology to design the higher-dimensional digital chaotic 
systems (HDDCS) in the domain of finite precision. The proposed 
system is based on the chaos generation strategy controlled by 
random sequences. It is proven to satisfy the Devaneys definition 
of chaos. Also, we calculate the Lyapunov exponents for HDDCS. 
The application of HDDCS in image encryption is demonstrated 
via FPGA platform. As each operation of HDDCS is executed in 
the same fixed precision, no quantization loss occurs. Therefore, 
it provides a perfect solution to the dynamical degradation of 
digital chaos. 

Index Terms —Chaotic encryption, Devaney’s chaos, dynamical 
degradation, FPGA implementation, high-dimensional digital 
chaotic system, random number generator. 


1. Introduction 

C HAOTIC system of simple structure can demonstrate 
complex dynamical properties in infinite mathematical 
world, such as sensitivity to initial conditions, topological 
transitivity and mixing, expansiveness, and decaying autocor¬ 
relation function ID. Such properties have subtle relation with 
some requirements of secure encryption system, especially 
sensitivity with respect to change of secret key O, O, 
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ll4l . So, designing chaos-based encryption schemes emerged 
as a new research direction to reinforce information secu¬ 
rity of data sent through the Internet O, IH, 171, HI, (91. 
However, any digital chaotic system implemented in finite- 
precision devices definitely degrade of various extents cni, 
nn . which may cause serious security flaws HJ), flTl . lfT4l . 
Multiplication, division, and other linear/nonlinear operations 
of real numbers may incur error by rounding, quantization 
and/or overfiow, which may further result in large differences 
between theoretical chaotic orbits and actual ones ca, ca. 
To counteract such degradation, many methods were proposed, 
such as perturbing chaotic states ifTTl . (TSl, perturbing control 
parameter CD, EqI, ED, or cascading two and more chaotic 
maps (22I, (231 , (24l . Recently, Deng and Hu showed that 
phase space of chaotic system must not be implemented in 
infinite precision, i.e., the definitions of chaos, like Devaney’s 
chaos, are compatible with domain of finite precision ini. 

Some works were done to disclose degradation character¬ 
istic of some simple digital domain chaotic maps, such as 
Logistic map (25l and piecewise linear chaotic map (26l. But, 
there is still lack of systematic theory to rigorously analyze 
and overcome degradation caused by finite precision effects. 
In 2010, Bahi et al proposed an 1-D integer domain chaotic 
system (IDCS) designated as Chaotic Iterations (Cl) system 
E3 The Cl system receives one or two random stream(s) 
as input and mixes them with chaotic iterations, where the 
involved mixture function is only bitwise OR operation. It 
was theoretically proven that the proposed systems satisfy 
the Devaney’s definition of chaos given in (28l . Since the 
systems run on finite sets of integer domains, the above 
problem caused by finite-precision disappears. As there is no 
any transformation from real numbers to binary sequences, Cl 
system can be considered as one of the most effective solutions 
for solving the degradation problem of digital chaotic systems. 

In (29l , a mark sequence was applied to avoid wasteful du¬ 
plication of values, leading to an obvious speed improvement. 
In (H, chaotic combination of two input entropic streams 
has solved flaws exhibited in the system designed in (29l . 
The chaos generation strategy was implemented in (SI, (321, 
( 33 II , through a sample-hold circuit and a decoder circuit 
so as to convert the uniform noise signal into a random 
sequence. The second chaos generation strategy named chaotic 
bitwise dynamical system (CBDS) was proposed in M- 
All the previous related works deal with 1-D chaotic maps, 
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the counterexamples in higher-dimensional domain and their 
dynamical properties are still unknown to date. 

In this paper, we study the general problem of constructing 
higher-dimensional digital chaotic system (HDDCS) in digital 
devices with finite precision, and establish general framework 
of composing HDDCS. We demonstrate that the associated 
state network of HDDCS is strongly connected, and then 
further prove that HDDCS satisfies the Devaney’s definition of 
chaos in the domain of finite precision. In addition, Lyapunov 
exponents for such a HDDCS are also given, and FPGA-based 
implementations are then provided to show the potential ap¬ 
plication of HDDCS in the digital world. The basic advantage 
and practical benefits of the proposed approach is that the 
chaos satisfying Devaney’s definition can be generated in the 
real digital devices by means of the chaos generation strategy 
controlled by random sequences. 

The remainder of this paper is organized as follows. The 
description of HDDCS is given in Sec. Jllj while proof of 
its chaos properties is provided in Sec. The Lyapunov 
exponents of HDDCS are calculated in Sec. W Section [V| 
presents application of HDDCS in image encryption with 
FPGA implementation. The last section concludes the paper. 

H. Design of HDDCS 
A. Description of HDDCS 

The general form of iterative equation for m-dimensional 
digital system with A'-bit fixed precision is 
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and P F Q = N. 

Set the general expression of m one-sided infinite random 
sequences as 


( 2 ) 


1 2 

n 

S= S S . . 

. S . 

1 2 

n 

b 

b 

II . 

. . U 

1 2 

n 

b 

b 

II 

. . V 


Likewise, each random number of these m sequences is 
expressed in binary form as 
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where n G Z+, V j G {P — 1, P — 2, ..., 0, —1, —2, ..., —Q}, 
s'j, u'j, and are the j-th bit of the binary form of 
and respectively. Note that the range of random numbers 
satisfy s'^^ G [0,2^ — 2“^] due to limited accuracy 
of A-bit representation. 

Through the iterative update mechanism controlled by ran¬ 
dom sequences, the general form of m-dimensional digital 
chaotic system (m-D DCS) is 
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denote the i-th component of the iterative function, where i = 
P-l,P-2,...,0,-l,-2,...,-Q. 

In Eq. we define chaos generation strategy controlled 
by random sequences for HDDCS as 
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if n" = 0, 


where j = P — 1, P — 2,..., 0, —1, —2,..., —Q. 

Then, Eq. © is further expressed as the following general 
form: 

'x^=ixy^ ■^) + {Fy)-s-), 


x^={xy^-TF) + iF2i-)-u^), 
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where the operators “(•)”. and “ + ” denote bitwise AND, 
bitwise NOT (negation), and bitwise OR, respectively. 
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Let us define £ as the set of points E of form 
((s, 14,..., i;), (xi, X 2 ,..., Xm)), where s, i4,..., are m in¬ 
dependent random sequences and xi, X 2 , . •Xm are A^-bit 
real numbers. Consider a metric space (f, d) and a continuous 
function Gp ^ 

Gpi^E^ G 14, . . . , u), {x \, ^2, • • • , ^m)) 

=((cr(s), (t(m), ..., cr(i;)), {xi,X 2 , 

Xm)), Hf 2 {i{u), {xi,X 2 ,..., Xm)), ■■■, {i{v), 

{Xi,X2,...,Xm))), 

( 6 ) 

where a{w) (w G {s,u, - ■ ■ , v}) shifts one integer in the one¬ 
sided infinite sequence w = {w^w'^ to the left, and 

'Hp^iis), (xi,X 2 ,...,Xm)) =((xi ■ iis)) + (Fi(-) • z(s))), 
Hp^(i{u), {Xi,X2, . . .,Xm)) =((X2 ’ i{u)) -h (-F2(-) ' *(«)))) 


MF^{i{v), {xi,X 2 , ■ ■ .,Xm)) = {{Xm ' i{v)) (-F’m(-) ' i{v)))- 

(7) 

As for function (j{w), one has 
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Fig. 1. The flowchart of the proposed higher-dimensional digital chaotic 
system. 
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Then, one can further get 

i{w) = 14 ;^, 

where w G {s,i4, ••• ,v}, k G and i{w) is equal to 
the overflow from the left shifting of the sequence w. In 
other words, some bits are randomly updated by Fi in each 
updating iteration, which are determined by 1 ( 5 ). Similarly, 
1 ( 14 ) determines the specific bits which are updated by F 2 , 
and i{v) determines the specific bits that are updated by Fm- 
The relation among m-dimensional digital domain system, m 
random sequences and m-dimensional digital chaotic system 
is shown in Fig. Note that true random sequences usually 
cannot be obtained via digital simulation 1^ , (SEl . But the 
input m random sequences can be generated by true random 
number generator (TRNG) based on digital circuit artifacts, 
such as metastable circuits (371, and digital oscillator rings 
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denote another one-sided infinite random sequences. Then, we 
introduce the distance d in the metric space (f, d) as 
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denote the k-th and (/c + l)-th iteration, respectively. 

With the above notations, HDDCS can be defined as 

E^^^ = Gp{E^), 

where /c = 0,1, 2,3,.... 
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Xf, Xf, X 2 , ^ 2 ,..., Xm^ Xm binary forms of real numbers 
with A^-bit finite precision, and 0 <d^ < pm{2^ - 2 ‘S)). 

B. Comparative study of RDCS, IDCS, CBDS, and HDDCS 

The proposed system HDDCS is briefiy compared with 
RDCS, IDCS in 02 and CBDS in 01, in Table ^ More 
detailed differences are summarized below: 

• The adopted implementation environment of RDCS is an 
infinite precision. However, IDCS, CBDS, and HDDCS 
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TABLE I 

Comparison of the digital chaotic systems, RDCS, IDCS, CBDS, and HDDCS. 


Name 

RDCS 

IDCS CBDS 
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Finite {N=P+Q , P bits integer, Q bits fraction) 

Dimension 

m-D 

1-D 

m-D 

External control inputs 

Null 

s = 

s" G {0,1,2,- 

s = •••s” ••• 

S — (%_i‘S'jv_2 ■■■‘^ 0)2 

x" e[0,2'"-l] 


S = ■■■s'" 

u = •••«"••• 

V = v'v^ ...y" ••• 

S = {,Sp_^Sp_2 '"Sq ■ S_^S_2 • • • S_Q {2 

u" = {Up_-^Up_ 2 '"K • w"ip 2 ■■■^- 2)2 

v" = iv;_yp_2 ■■■Vo- ^>"-2 • • • v"e >2 

g[0,2^-2“2] 

Iterative equations 

A ■ 

II II II 

rF(p-'). if 7=p 
x" =\ 

\xC if j * s'" 

j = N-\,N-2,---,0 

if =1 

x" =< 

xC if s" = 0 

j = N-l,N-2,---,0 


fF,0. if^;=i 

0=^1 (•).»= / 

' [xi J if Sj = 0 

fp(-). ifM"=l 

^ {X 2 J if Uj = 0 

\f (•). if v" =1 

<y=P(-L= „ 

' Wj if U = 0 

/ = P-l,P-2,---,O,-l,-2,---,-0 


are implemented in finite precision, they are suitable for 
digital computers or other digital devices. 

• RDCS present the finite integer part and infinite fractional 
part of a real number as follows: 
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But IDCS and CBDS process data in integer domain as 
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Furthermore, in HDDCS, the range of the input data are 
extended to digital domain of finite integer part and finite 
fractional part as 
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where N = P ^ Q. 

• Both IDCS and CBDS deal with 1-D chaotic systems. 
In contrast, HDDCS extends the system to any finite 
dimension. 

. Compared to RDCS, IDCS, CBDS and HDDCS need 
external inputs. 

• The main features of discrete-time RDCS is that all the 
bits in Xn-i will be updated by iteration function F 
at each iteration operation. Likewise, all the bits in Xn 
are updated by iteration function F at each iteration 


operation. But only one bit of Xn in IDCS is updated by 
iteration function F at each iteration operation. CBDS 
uses multiple random bitwise operations instead of only 
one in IDCS. HDDCS are similar to CBDS but can work 
in higher-dimension. 

C. Network analysis of the state space of HDDCS 

Given a digital chaotic system, a state and its interval is 
mapped to another one. Considering each state or interval 
as a node (vertice), and the mapping relation as a directed 
edge (link), state network of the chaotic system can be build 
up. As shown in 13^ . iQi, the associated state network 
can demonstrate some dynamical properties of digital chaotic 
systems that can not be observed by the previous analytic 
methods. For the directed graph of Gp in HDDCS, all the 
possible combinations of (xi, X 2 ,..., Xm) are the nodes, and 
there is a directed edge from node (f i, ^ 2 ,..., Xm) to another 
node {xi,X2 ,... ,Xm) if 

{Gpifs^U^ . . . ,i)), (^1, X2-) • • • 5 XjYi)y)x\,X2^---iXm ~ 

{X \, X2 •)••••) XjYif 

We found that the state network of HDDCS must be strongly 
connected, namely every node is reachable from any another 
one. 

In the following, we use 2D-DCS with N = 2 {P = 2, Q = 
0) to illustrate the property on connectivity. For example, con¬ 
sider a 2-D digital system uncontrolled by random sequences, 

as ^ _ 

(x^ = = x--\ 

\y^ = F2{x^-\y^-^) 
where 0 denotes the bitwise XOR. 























YU et al: 


5 


TABLE II 

State transition for system 



(x",2/") 


(x",y") 
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(0.1) 

(1,2) 

(2,0) 

(3,2) 

(0.2) 

(1,3) 

(2,1) 

(3,3) 

(0.3) 


The concrete state transition is shown in Table [TIJ and the 
corresponding state transition diagram is shown in Fig. 




Fig. 2. The disconnected network for Eq. 

According to Eq. (0, Eq. 0, and Eq. ([8]), one can obtain 
the iterative form of GF{E)x,y as 

where s = s^s‘^s^ ... and u = ... are two random 

sequences generated by TRNG. 

According to Eq. ([^, the concrete state transition of 
GF{E)x,y is shown in Table H and the corresponding state 
network is shown in Fig. , which is a strongly connected 
network. In the next section, we will prove that HDDCS 
satisfies the Devaney’s definition of chaos if the state network 
of its Gi? is strongly connected. 



Fig. 3. The strongly connected network for of system ([8|. 


Proof: If w'^ = w'^ for i = 1, 2,..., n , then 
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Due to the definition of the proposed distance: for any m ^ n, 
if = w'^, then d{w,w) < ■ 

The lemma can let us quickly determine whether the two 
sequences are close to each other. From intuitive observation, 
we can assure two sequences are close to each other as long 
as they have a considerable number of consistent foregoing 
entries. 


III. Chaotic performance of HDDCS 

In this section, we prove that the map of HDDCS satisfies 
Devaney’s definition of chaos, i.e., its periodic points are dense 
in its definitional domain and it is transitive ED. 

A. Dense periodic points of HDDCS 

To prove the existence of a dense set of periodic points 
and the transitivity property of the discrete dynamical system 
defined previously, the following lemma is firstly proven. 


Theorem 1. The periodic points of HDDCS are dense in the 
metric space {S^d). 

Proof: For any given 5 G (0,1) and a periodic 
point ((s, ..., X 2 ,..., G S, let us try to 

prove that the dense periodic points in {S,d) can always 
be found in the neighborhood of distance 5 of any point 
{{s,U, {xi,X2,.. .,Xm)) e ^ as 

d{{{s,U, ■ ■ ■ ,v),{Xi,X2, . . . ,Xm)), 

{{S,U, . . . ,V),{XI,X 2 , . . . ,Xm))) < S- 


Lemma 1. Let w G {s, it, • • ■ , v}, w = ... it;" ... 

and w = w^uP"ufi ...«)" ..., the metric distance d satisfies 
that if ID® = It)® for i = 1,2,3 ,... n, then 

d{w,w) < 

where ,w^ G [0, 2^ — 2~'^] for k G Z+. 


Without loss of generality, we assume that the general form 

of ((s,M, {xi,X2, ■ . ■,Xm))) is 

{{S,U, ... ,V), {XI,X 2 , ... ,Xm))) 

= . . . S" . . .), it" .. .), 

. . . , {v^v^ ...v'‘°...v^.. .)), {xi,X2, . . .,Xm)) £ £■ 
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TABLE III 

State Transition of Gp for system ^ in HDDCS. 


/ n -1 

(S",-!!") 

\x : y ) 

(0,0) (0,1) (0,2) (0,3) (1,0) (1,1) (1,2) (1,3) (2,0) (2,1) (2,2) (2,3) (3,0) (3,1) (3,2) (3,3) 

(0,0) 

(0,0) (0,1) (0,2) (0,3) (1,0) (1,1) (1,2) (1,3) (2,0) (2,1) (2,2) (2,3) (3,0) (3,1) (3,2) (3,3) 

(0.1) 

(0,1) (0,0) (0,3) (0,2) (1,1) (1,0) (1,3) (1,2) (2,1) (2,0) (2,3) (2,2) (3,1) (3,0) (3,3) (3,2) 

(0.2) 

(0,2) (0,3) (0,0) (0,1) (1,2) (1,3) (1,0) (1,1) (2,2) (2,3) (2,0) (2,1) (3,2) (3,3) (3,0) (3,1) 

(0.3) 

(0,3) (0,2) (0,1) (0,0) (1,3) (1,2) (1,1) (1,0) (2,3) (2,2) (2,1) (2,0) (3,3) (3,2) (3,1) (3,0) 

(1.0) 

6,0) 6,0) 6,2) 6,2) 6,0) 6.0) 6.2) 6.2) 6.0) 6.0) 6.2) 6.2) 6.0) 6.0) 6,2) 6,2) 

(1.1) 

6,1) 6,1) 6,3) 6,3) 6,1) 6.1) 6.3) 6.3) 6.6 6.6 6.3) 6.3) 6. i) 6. i) 6,3) 6,3) 

(1.2) 

6,2) 6,2) 6,0) 6,0) 6,2) 6.2) 6.0) 6.0) 6.2) 6.2) 6.0) 6.0) 6.2) 6.2) 6,0) 6,0) 

(1.3) 

6,3) 6,3) 6,6 6,6 6,3) 6.3) 6.6 6.6 6.3) 6.3) 6.6 6. i) 6.3) 6.3) 6,1) 6,1) 

(2.0) 

6,0) 6,1) 6,0) 6,6 6,0) 6.6 6.0) 6.6 6.0) 6.6 6.0) 6. i) 6.0) 6.1) 6,0) 6,1) 

(2.1) 

6,1) 6,0) 6,6 6,0) 6,6 6.0) 6.6 6.0) 6.6 6.0) 6.6 6.0) 6.1) 6.0) 6,1) 6,0) 

(2.2) 

6,2) 6,3) 6,2) 6,3) 6,2) 6.3) 6.2) 6.3) 6.2) 6.3) 6.2) 6.3) 6.2) 6.3) 6,2) 6,3) 

(2.3) 

6,3) 6,2) 6,3) 6,2) 6,3) 6.2) 6.3) 6.2) 6.3) 6.2) 6.3) 6.2) 6.3) 6.2) 6,3) 6,2) 

(3.0) 

6,0) 6,0) 6,0) 6,0) 6,0) 6.0) 6.0) 6.0) 6.0) 6.0) 6.0) 6.0) 6.0) 6.0) 6,0) 6,0) 

(3.1) 

(3,1) (3,1) (3,1) (3,1) (2,1) (2,1) (2,1) (2,1) (1,1) (1,1) (1,1) (1,1) (0,1) (0,1) (0,1) (0,1) 

(3.2) 

(3,2) (3,2) (3,2) (3,2) (2,2) (2,2) (2,2) (2,2) (1,2) (1,2) (1,2) (1,2) (0,2) (0,2) (0,2) (0,2) 

(3.3) 

6,3) 6,3) 6,3) 6,3) 6,3) 6.3) 6.3) 6.3) 6.3) 6.3) 6.3) 6.3) 6.3) 6.3) 6,3) 6,3) 


Given 5 < 2 “^, if point (xi, X2,..., and 

(xi,X2,..., Xm) do not coincide in the m-dimensional space, 
we can obtain (xi, ^2,..., Xm) 7^ (^1, ^2, • • •, ^m) such that 

• • • 5 ^m)i (^l5 ^2', • • • 5 ^m) 

= U (*1 - *1)^ + (*2 - X 2 y + ... + {Xm- XmY 
> 2~^. 

Then, one has 

d(((s,U, (xi,X2, . . ■,Xm)), 

((§,u,...,v),(xi,X2,...,X^))) > e. 

So, to satisfy 

d(((s, u,..., v), (xi, £2,..., Xm)), 

((s,U,...,v),(xi,X2,...,Xm))) < £, 
we must first set xi = xi, £2 = £2,, Xm = In case 

d(((s, U,...,v),{xi,X2,..., Xm)), 

{{s,U, . . . ,v),{xi,X2, . . . ,Xm))) <£, 

we should consider to prove that 
((5, fi,..., i)), (xi, ^2,..., is a periodic point. If 

the first elements of s and s are the same, then 
ds{s,s) < < 5. Referring to Lemmasimilar results 

can be obtained as ko elements of , i) and fi,..., are 
the same and 


The above equation shows that HDDCS starts from 
(xi, X 2 ,..., then returns back to it after ko iterations. 
This means that a periodic point (f 1 , ^ 2 ,..., is found 

out: 

{{s,u, ...,v), {xi,£ 2 , ■ ■ .,Xm)) 

= ...), {u^u^... u^°v}u^ ... u’^°. . 

.)), {xi,X 2 , .. ■,Xm)), 

which satisfies 


G'p{{s, U,...,v), {xi,X2, Xm)) 

= {{s,u, . . . ,v),{xi,X2, . . . ,Xm)), 


d{{{s,u, ...,v), {xi,X2 , . . . , Xrn))’) 

{{s,U,...,v),{xi,X2,...,Xm))) 

= ds{s, s) + du{u, u) + ... dy{v, v) 

^2? • • • 5 (^1? ^2? • • • 5 ^m)) 

= ds{s, s) + du{u, u) ^ dy{v, v) 

< s. 


If after ko-th iteration, one get 


{G'p{(s,u, 


and 


{xi-, X2i • • • 5 ^m)))cCi,CC2v?3:m 

^ (xi,X 2 , . . .,Xm) 


du{u, u) <2 < s,... ,dy{v,v) <2 < 5. 

So, V 5 < 1, an integer ko satisfying the relation ds{s,s) + 
du{u^u) + ... + dy{v,v) < m X 2“^^° < 5 can always be 
found. For instance, to make m x 2“^^° < 5 hold, the value 
of ko can be set as 

ko = L(log 2 rn - log^ s)/N\ + 1. 

After ko-th iteration, one has 


, X2, • • • 5 Xm)))xi,X2,---,Xm 

(^l5 ^25 • • • : ^m)' 

Because is strongly connected, there is at least one path 
from the state {x[,X2, ..., x'^) to the state (f 1 , ^ 2 ,..., Xm), 
after another iteration of io times, where io is equal to the 
number of edges in the connected path between the state 
^ 2 ,..., and the state (xi, X 2 ,..., Xm)- By making 
the equation 


, X 2 , • • • 5 Xm)))xi,X 2 ,---,Xm 

— (^1 5 X 2 ^m)' 


...,V), {XI,X2, . . . ,Xm)))xuX2,...,Xm 

— {x\ , X2 5 • • • 5 ^m) 
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hold, a periodic point is found by checking 


((s,n, 

v), {xi,X2, . . ■,Xm)) 


gko ^ko + 1 ^ko +2 ^/co+io^l^2 ^ko 

gko+lgko +2 ^ ^ ^ gko+io _ ^ (^1^2 ^ ^ ^ ^t^o^feo+l 


.. ... 

y^ko+io 

+ . . . 

yko-\-io^ 

,1^2 ykoyko-\-lyko -\-2 yko-\-io 

{XI,X 2 , 

• • • 5 Xm))-, 

which satisfies 

d{{{s, u. 

,...,v),{xi,X2,...,Xm))), 

{{S,U, . . . ,V),{XI,X 2 , . . . ,Xm))) 


= ds{s, s) + du{u,u) + ... + dy{v,v) 

~^dx{^X\ 5 X2 5 • • • : Xra)’) {x\ , X2 5 • • • 5 ^m) ) 

= ds{s, s) + du{u,u) + ... + dy{v,v) 

< e. 

In summary, the periodic points of are dense in metric 
space as shown in Fig. ■ 



Fig. 4. The diagram of the periodic points of HDDCS. 

B. Transitive property of HDDCS 

Theorem 2. Function Gp is topological transitive in the 
metric space (f, d). 

Proof: The so-called topological transitivity of function 
Gf in metric space means that there always exists 

no > 0 satisfying G^{U') H U" 0 for any nonempty 
open sets U' and U". Let {{s'^u'^... ,7;'), {x'i^x' 2 ^ ... ^x'^)) G 
U' C E denote the center of U', which can be represented as 



),(Xi,X2,..., 

X'm)) 

II 

... s'"“... s'" 

...), 

(n'U'2 . 

.. ... m'" . 



.. v'^° ... n'" . 

••)), 

(xi,X2,, 

. . . , Xm)' 



The center of U", {{s", 

to 

also be presented as 



((s",n",..., 

v"),ix'l, 

X 2 , ■ ■ ■ , Xm 

II 

... 5 ""“. 

..s""...), 


. ... 

n"".. 

(u"U"2., 

, . .. 


-r" t" 

5 *^2 5 • • 

x" i 

. , ^mJi 



and a point in U' is denoted as 
(( 5 , n,..., i)), (xi, X 2 ,..., Xm) eu' C£. 

If the sphere radius of U', r', is less than 2~^ and points 
{x'i^x' 2 ^..., x'^) and (xi, X 2 ,..., Xm) do not coincide in the 
m-dimensional space, one has 

{{s,U,...,v),{xi,X 2 ,...,Xm) i U' 

as 

dx{{xi,X2, Xm), {x'i,x'2, ■ . ■ ,x'm)) 

= \J (*1 - + (*2 - 4)^ + ... + {Xm- x'mY 

> 2“^ > r', 

and (xi, X 2 ,..., Xm) 7 ^ (^ 1 : ^2 5 • • • 5 ^m)- to satisfy 

((s,M, . . . , v), (^1,52, • • ■,Xm) e U', 

we must set 

(xi, X2’, • • • : Xm) (^1: ^2’ • * • 5 ^m) 

to obtain 

dx{{xi,X 2 , . . .,Xm), {x[,X 2 , X^)) = 0. 

If the first ko elements of s' and s are the same, then 
ds{s', s) < From Lemmathe similar results can be 

obtained when the first ko elements of u',... ,v' and fi,..., i) 
are the same and 

du{u', u) < < 5,..., dy{v', v) < 2-^^° < 

So, V r' < 1, an integer ko satisfying ds{s', s) + du{u',u) + 
... + dy{v', v) < m X 2“^^° < r' can always be found. For 
instance, to satisfy m X 2-^''° < r', the value of ko can set 
as 

ko = L(log 2 m - log 2 r')/N\ + 1. 

If after the ko-th iteration, equality 

, X2, . . . , Xm)))xi,X 2 ,---,Xm. 

(Xi, X 2 , • • • , Xm) 

exists, then no = fco is found, and 

{{s,U, . ■ ■ ,v), {xi, X 2 , ■ . ■ ,Xm)) 

= ,,, s"” ...), (n'U'2 ... 

... n"” (n'h'2 ... _ _ 

P^...)),{x[,x' 2 ,...,x'm))eU', 
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which satisfies 

Xm)) 

G G^p°iu')nu" 

So, 


hold. 

If after the ko-th iteration, inequality 

{Gp° {{S, U,...,v),{xi,X2,...,Xm))) 

Xi,X2,---,Xm 

^ {x'{,x' 2 ,...,x'^), 

holds, set {G'P ({s,u,... ,v), (xi, X 2 ,... ,Xm)))xi,x 2 ,...,x^ = 
(x'/', ^ 2 Because is strongly connected, there 
is at least one path from the state {x'{',X 2 ',... to the 
state (x'/, X 2 ,..., x'^) after other io iterations. As 

hold, riQ = ko io is found, and 
{{s,u, {xi,X2, . . .,Xm)) 



gfko ^/co + l^/co+2 

gko+io ^nl ^112 ^ 


^0+2 

1X"U"2... 


^ko + l 

...V 

,feo+jo^"l^"2 _ _ _ ^//n 

•••)), 

(x{,x'2,... 




e U, 

which satisfies 

Gp°{{s,U, {xi,X 2 , . . .,Xm)) 

e G^p°{U')nu''. 


So, one has 

In summary, G^ is transitive in the metric space {S^d). ■ 

If a dynamical system on a metric space is transitive and 
has dense periodic points, then it has sensitive dependence on 
initial conditions ED- In other words, if the state network of 
HDDCS is strongly connected, we can prove that it is chaotic 
in the sense of Devaney’s definition of chaos. 


IV. Lyapunov exponents of a class of HDDCS 


A. The general expression of equivalent decimal for Gp 

Set the binary form of m-dimensional array of N-bit 
integers as 


'xi = Xi^N-1Xi,N-2 • • • ^ 1 , 0 , 

X2 = ^2,Ar-1^2,Ar-2 • • • ^2,0, 

^m,N—l^m,N—2 • • • ^m,0: 


where Xij G {0,1} (i = 1, 2,..., m; j = N—1, N—2, ..., 0). 

The general form of corresponding decimal integer is ob¬ 
tained as 

/ AT-l 

k=0 


X2 


N-1 


k=0 


( 10 ) 


AT-l 


k=0 


Random number of m sequences is expressed in correspond¬ 
ing decimal fraction as 


' +00 

s= E(s'' • 2“^''), 

k=l 
+ 00 

u= Eu'' •2“^'=), 

k=l 


( 11 ) 


+ 00 

v= 

k=l 


where 2~^^{k = 1,2,...) denote weights. 

According to Eq. ( p^ and ( pTj ) added together, the general 
form of the corresponding decimal number is obtained as 


N—l +00 

2/1 = E(^i.fe-2")+E("''-2-^")’ 

k=0 k=l 

N—l +00 

2/2= E(^2,;c-2'=) + EU''-2-'^"), 

\ k=0 k=l 

N—l +00 

ym= E(^-.'i-2") + E(^''-2”'^")- 

^ k=0 k=l 

According to the chaos generation strategy controlled by 
random sequences, the general form of m-Dimensional digital 
discrete-time iterative equation can be presented as 


'fli(X) =(x,-s0 + {FX)-s^), 

g2{X) =(a;2-U) + (F2 (-)-U), 


.gm{X) = {Xm ■ vO + {Fm{-) ' V^), 


In this section, the Lyapunov exponents of HDDCS with 
N = P {Q = {)) are estimated. 




(12) 
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where X = (Xi, X2,..., X^) and 
^Fl(-) =Fi{xuX2, 

F2{’) =F2(xi,X2, 

< 

< Fmi,'^ —Ffn (^ 1 1 X2 : • • • 5 ^m) • 

Separately shifts one value in each one-sided infinite se¬ 
quence (s = 8^3“^ . . . U = U^V? . . .U^ . . V = 

v^v‘^ .. .v'^ ..the first value turn into v‘^ indi¬ 

vidually, the corresponding weight is 2 “^. So, general form 
of the corresponding decimal fraction after shifting one value 
in every one-sided infinite sequences is 



+00 


g,iS) 

k=2 

+00 

g2{U) 


^-Nk\ 


k=2 


+00 


(13) 




k=2 


Adding Eq. and (HI* together, one can obtain the 
general form of the corresponding decimal number. 


9i{yi,y2,---,ym) 

= (xi • + (Fi(-) ■ fii) + 2^ 

k=2 

52 ( 2 / 1 , 2 / 2 , •••,2/m) 

+ 00 

= {X2 ■ + {F2{-) ■ U) + 2^ • 2-^'=), 


+00 


k=2 


Qmiyi-) y2’s • • • 5 Vm) 

+00 

= {Xm ■ + (F„(-) ■ V^) + 2"^ • 2“'^''). 

, k=2 

after randomly updating multiple random bits and shifting one 
value in every one-sided infinite sequences. 

B. The mathematical expression for ^3i{yi,y2,---,ym) 

In the interval [^, (n G [0,2^^ ~ 1]). the part of 
decimal integer is not changed that AXi = AX 2 = • • • = 
AX^ = 0. Furthermore, the first decimals of the m sequences 
are the same, so 

+00 

At/i = AXi + A5 = A5 = ^(As'' • 2"^*), 

k=2 
+ 00 

Ay2 = AX 2 +AU = AU = ^(Au'' • 2“^''), 

\ k=2 

+00 

Ay„ = AA„ + AU = AU = • 2-^'=). 

, k=2 


From the definition of partial derivative, one has 
dgi{yi,y 2 ,---,ym) 


dyi 

lim 

Ayi^O 


lim 


gi(5i +Ayi,y 2 ,--- ,ym) - 5i(5i, 52, • • • ,5m) 
^ Ayi 

gi(A) + 2^Etr2((g'^ + A.")-2-^'=) 

ES(a. 


= lim 


2^ 




N 


= 2 

Similarly, one can get 

dg2{yi,y2,--- ,ym) dgs{yi,y2,--- ,ym) 


dy2 


dyz 

dg m ( 51 ,52,--- ,5m) 


dy„ 


= 2 


N 


On the other side, one has 


951 ( 51 , 52 ,-- 

5 ym) 

952 


952 ( 51 , 52 ,-- 

•) ym) 

951 


95m(5l,52,-- 

•) ym) 


dyi 


%i(5i, 52,-- - ,5m) ^ _ _ 

dys 

dgi{yi,y 2 ,--- ,ym) ^ „ 
dym 

^ 52 ( 51 , 52 ,-•• ,5m) ^ ^ ^ 

dyz 

^ 52 ( 51 , 52 ,-•• ,5m) „ 

dym 

m (51,52,-•• ,5m) ^ ^ 

dy2 

m ( 51 , 52 ,-•• ,5m) ^ g 

%m-l 


So, the corresponding Jacobian matrix, diagonal matrix 
diag(2^, 2^, • • • , 2^), is obtained. 


C. Estimating the Lyapunov exponents 

Let • ^n) denote the k-\h characteristic value of 

matrix (<I>^ • the Lyapunov exponent of the specific 
HDDCS can be estimated as 


KVk) = 


lim — 

n^+oo 2n 

lim — 

n^+oo 2n 


ln|Al,(C-^n)l 

ln((2^)2^) 


= Xln2, 


where and /o = 1,..., m 14^ . 


V. FPGA-based real-time application of 3D-DCS 

In this section, we first give the design and implementation 
of an FPGA-based generator for 3D-DCS. Then, a RGB color 
image encryption method based on 3D-DCS is presented. We 
split the RGB color image into its three R, G, B components. 

























10 


IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS-I: REGULAR PAPERS, VOL. 63, NO. 3, MARCH 2016 


and use the 3D-DCS to scramble the pixel values of the three 
components R, G, B. Finally, a systematic methodology for 
FPGA platform-based implementation of the above method is 
proposed. 

A. Design of 3D-DCS in FPGA 

To utilize the power of FPGA, the computation needs to 
be divided in several independent blocks of threads that can 
be executed simultaneously. The performance on FPGA is 
directly related to the number of threads and that of logistical 
elements used during processing, its performances decrease 
when number of branching instructions (moves like while, if, 
etc.) increases. Following these rules, it is possible to build a 
Verilog-HDL program of the 3D-DCS algorithm 1431 . 

A concrete example is provided to illustrate the 3D-DCS 
algorithm. Here, we use 3D-DCS with N = 32 (P = 32, Q = 
0). For example, a 3-D digital system controlled by random 
sequences is considered: 

^x'^ = 0 (1 « (mod( 2 :’^“^, 32))), 

< = yn-l 0 (1 « (mod32))), (14) 

=^^0(1 « (mod(^^-\32))). 

where mod(m, n) = m — n • [^J, gives the largest 
integer less than or equal to and « represents the left 
bit shift. According to Eq. (©,Eq. (|7]) and Eq. ( p^ , one can 
obtain the iterative form of GF{E)x,y,z as 

= x^-Pp^0((^^0(l « (mod(z^-\32))))-5^), 

' y"" « (mod(x^“\32))))-'u^), 

« {mod{y^-\32))))-v^). 

(15) 

where s = s^s‘^s^ ..., u = ... and v = v^v‘^v^ ... 

are three random sequences. This 3D-DCS may utilize any 
reasonable random sequence as s^u^v. For demonstration 
purposes, ISAAC in ll44]| is adopted here. It was found that 
the three channels of the outputs of 3D-DCS all can pass 
the NIST randomness test suite given in ia. Meanwhile, the 
correlations among the three channels of the outputs and their 
auto-correlation strengths are all very low. 

Figure depicts the circuit structure of 3D-DCS. First 
of all, according to (381, three different types of oscilla¬ 
tor rings TRNGs named as Oscillator_Ringsl:instl, Oscilla¬ 
tor _Rings2:inst2 and Oscillator_Rings3:inst3 blocks shown in 
Fig. are applied as the external control inputs s, u, and v. 
Then, the block 3D-DCS_jyrocessing:inst4 is constructed by 
using Eq. At last, the states of oscillator rings TRNGs 
and 3D-DCS_j)rocessing:inst4 are updated in the feedback 
block controlled by clock signal, and it executes in clock 
positive edge. The frequency of input clock decides the speed 
of the 3D-DCS processing. In our experiments, the input clock 
is set as 50MHz. Then, ModelSim Altera is used to obtain 
= outl, y'^ = out2 and z'^ = out3, which are used for 
RGB image encryption and decryption processing. 

B. Design of the FPGA-based hardware system for image 
encryption and decryption 

Figure shows the block diagram for FPGA-based hardware 
system of image encryption by 3D-DCS, where the function 



Fig. 5. Block Diagram of 3D-DCS in FPGA. 


of the FPGA hardware part is implemented by Verilog-HDL 
program. 



Control Signal 


Fig. 6. Block Diagram of FPGA-Based application on image encryption. 

As shown in Fig. the hardware system at the transmitter 
side consists of five parts: picture RAM, 3D-DCS, VGA 
display controller, monitor 1, and RGB encryption module. 
The picture is encrypted at the transmitter side and then 
transmitted through the public channel to the receiver. The 
receiver, for its part, contains four components: monitor 2, 
VGA display controller, RGB decryption module, 3D-DCS. 

The corresponding hardware implementation platform is 
shown in Fig. |7j two same models of Altera DE2 FPGA 
development board were used. The working principle of the 
hardware system is expressed as follows. In Fig. [TJa), the 
picture is previously stored in the Altera DE2’s RAM (Random 
Access Memory), namely Picture RAM in Fig. When the 
system is turned on, the picture is sent in two ways: one is 
delivered to VGA display controller, it collects the RGB (Red, 
Green, and Blue) information of every pixel from the picture, 
then transmits them to monitor for display; the other way is 
sent to RGB encryption module, the RGB elements of each 
pixel are encrypted by operating XOR with the output stream 
of 3D-DCS, respectively. This needs three different types of 
oscillator rings as external input sources. In addition, the 
control signal of states for three oscillator rings at the sender 
are also delivered through the public channel to synchronize 
the ones of 3D-DCS at the receiver. 

In Fig. I^b), the received encrypted RGB objects are de¬ 
crypted. Since the states of the 3D-DCS synchronize with the 
ones at the transmitter side, the RGB values of the picture 
could be recovered by executing XOR separately with the 
states at three dimensions of 3D-DCS. The RGB pixels are 
finally sent to VGA controller and displayed on the monitor. 
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C) d) 

Fig. 7. FPGA-based implementation results for chaos-based secure image 
communications: a) Two Altera DE2 FPGA development boards; b) Hardware 
platform; c) Case of mismatched parameter; d) Case of matched parameters. 

C FPGA-based implementation result for image encryption 
and decryption 

In the experiment, an image of 640 x 480 resolution is 
previously stored in Altera DE2 board’s RAM. The image 
encryption system based on 3D-DCS is implemented with 
this image. As shown in Fig. only when the states of 
3D-DCS between the transmitter and receiver are exactly 
matched, then the original image can be correctly decrypted, 
otherwise it cannot be recovered. The experimental FPGA- 
based implementation results are shown in Figs. [TJc) and (d). 

VI. Conclusion 

To solve the degradation of chaotic dynamical properties 
caused by limitation of finite-precision presentation and quan¬ 
tization, this paper developed a novel higher dimensional 
digital chaotic system (HDDCS), utilizing the chaos generation 
strategy controlled by random sequences. It is proved that 
HDDCS satisfies the Devaneys definition of chaos. Finally, 
FPGA-based implementation of image cryptography algorithm 
with 3D-DCS is detailed, confirming the feasibility and ap¬ 
plicability. This work provides perfect solution to dynamical 
degradation of digital chaos and may further promote appli¬ 
cation of chaos in digital world. 
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